Steering effects on growth instability during step-flow growth of Cu on Cu(l,l,17) 
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Kinetic Monte Carlo simulation in conjunction with molecular dynamics simulation is utilized to 
study the effect of the steered deposition on the growth of Cu on Cu(l,l,17). It is found that the 
deposition flux becomes inhomogeneous in step train direction and the inhomogeneity depends on 
the deposition angle, when the deposition is made along that direction. Steering effect is found to 
always increase the growth instability, with respect to the case of homogeneous deposition. Further, 
the growth instability depends on the deposition angle and direction, showing minimum at a certain 
deposition angle off-normal to (001) terrace, and shows a strong correlation with the inhomogeneous 
deposition flux. The increase of the growth instability is ascribed to the strengthened step Erlich 
Schwoebel barrier effects that is caused by the enhanced deposition flux near descending step edge 
due to the steering effect. 

PACS numbers: PACS numbers: 68.35.-p, 68.37.-d 



I. INTRODUCTION 

In the growth of thin fihus on a vicinal surface of high 
areal step density, there is a net current of deposit par- 
ticles towards the ascending step edge due to the Er- 
lich Schwoebel barrier at descending step edge. Such 
transport of deosit atoms to the ascending step increases 
the possibility of step flow growth and makes the growth 
of thin films on a vicinal surface more stable than that 
on a sigular surface ji*^ Moreover, such asymmetric flow 
of deposit atoms provides the possibility of forming a 
structure along the step edge and has been a subject 
of numerous studies for the growth of one-dimensional 
systems]^ Even in the thin film growth on a vicinal sur- 
face, however, develops the meandering instability along 
step-edge j^i^i^ Possible sources for this instability^ have 
been suggested as the asymmetric adatom diffusion due 
to the step Erlich-Schwoebel barrier^-^ and kinetically 
limited diffusion of the adatoms due to the kink Ehrlich- 
Schwoebel barrier 

In addition to the kinetic effects mentioned above, the 
deposition process, one of the ignored dynamic processes, 
has been recently found to affect the thin film growth^. 
That is, the interaction between a deposit atom and the 
atomic structure on the surface modifies the trajectory 
of the deposit atom, called steering effect, and causes 
the inhomogeneous distribution of adatoms affecting the 
growth of thin filmsiSaiS, Adjacent to the edge of islands or 
steps, the steering effect is conspicuous due to rapid vari- 
ation of the interaction potential, and thus is expected 
to be more influential in a deposition on a vicinal surface 
having high areal step density than in that on a singular 
surface. 

The purpose of the present study is to explore the role 
of steered deposition on the thin film growth on a vicinal 
surface, which has been ignored in most of the previous 
simulation or theoretical studies (see Ref. 2 for a re- 
view). Specifically, we calculate the deposition fiux dis- 



tribution on the vicinal surface, varying the deposition 
angle and direction. We also search for any possibility 
to overcome such kinetic growth instability by adjusting 
the dynamic variables involved in the deposition process. 
We have chosen to study the growth of Cu on Cu(l,l,17), 
because Cu(l,l,17) shows no surface reconstruction and 
has been a subject of many experimental and theoret- 
ical studiegSiii allowing us to compare our results with 
preexisting ones. Present study utilizes a computer simu- 
lation combining a molecular dynamics (MD) simulation 
for the dynamics of deposit atoms with a kinetic Monte 
Carlo (KMC) simulation for the growth of adatoms on 
the surface. 

We find that the steering-induced enhancement in the 
deposition fiux near descending step edge is a critical 
factor affecting the growth instability on vicinal surface. 
The inhomogeneity of deposition flux depends on de- 
position angle, and a deposition angle which gives the 
minimum growth instability is found. Nevertheless, the 
steering effect always increases the growth instabililty re- 
gardless of the deposition angle, with respect to the case 
where steering effect is neglected. 



II. SIMULATION METHOD 

KMC simulation is adopted to simulate the whole pro- 
cess of thin film growth. Contrary to conventional KMC 
scheme, to simulate the trajectories of depositing atoms in 
detail, we incorporate MD into KMC simulations, where 
MD is employed whenever a deposition event occurs in 
the KMGia. 

In the MD simulation, a Lennard-Joncs potential 
U{r) ~ 4D[{a/r)^'^ — (u/r)^] is used for the pair inter- 
action between a deposit atom and an atom on surface, 
with D = 0.4093 eV and a = 2.3381. These values of 
D and a are adopted from Dijken et.alMii^. The ini- 
tial kinetic energy of the deposit atom is set to 0.15 eV, 
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TABLE I: Difltusion barriers and parameters used in KMC. 
Same notation is used for each diffusion process as in Fig. 1. 



diffusion type 


diffusion barrier 


El 


0.42 eV 


E2 


0.38 eV 


E3 


0.51 eV 


E4 


0.68 eV 


E5 


0.59 eV 


E6 


0.18 eV 


ES 


El+O.l eV 


jump frequency (fo) 


3.6 X 10'^ 


deposition rate (Fq) 


0.003 ML/s 



corresponding to the melting temperature of Cu. The 
Newton's equation of motion is solved using Verlet algo- 
rithm. Atom is approached to the substrate by MD, and 
then positioned to the nearest four-fold hollow site from 
the terminal position. The transient mobility is not in- 
cluded in the present study. That is, the deposit atoms 
are assumed to be in equilibrium with the substrate right 
after the deposition. 

In the KMC simulation, a lattice gas model is adopted 
which allows jump diffusion for adatom motion on fee lat- 
tice. The possibility of each jump diffusion is calculated 
from the corresponding hopping rate, ly = voexp^^^ , 
with attempt frequency, lyQ = 3.6 x The def- 

initions of the most relevant diffusion processes in the 
present simulation are illustrated in Fig. 1. In Table I, 
listed are the values of the diffusion barriers, Ei, that are 
adopted from the values used by Koponer&ii in a growth 
simulation of Cu on Cu(l,l,17) and those obtained by 
Furmann^'^ from a simulation study for thin film growth 
on Cu(OOl). 

Cu(l,l,17) surface has a (001) terrace of 8.5-atomic 
width between two steps of an atomic height in [-1,1,0] 
direction. In the following, x-axis is along the step edge 
as shown in Fig. 1, and y-axis is along the step train 
direction. The simulation box has 12 terraces with step 
edge length of 800 oq, where uq is the surface lattice con- 
stant of Cu(OOl), 2.55 A. Periodic boundary conditions 
are adopted in both x and y directions. 





FIG. 1: Illustration of some diffusion processes taken into 
account in the present simulation. 



FIG. 2: Trajectory of deposit atoms and normalized deposi- 
tion flux, (a) Trajectory of deposit atoms at deposition angles 
of (a) —70° and (b)+70°. Steered deposition fluxes at depo- 
sition angles of (c) -70° , (d)+70°, (e) -35° , (f)+35°, (g) 
0°, and (h) deposition without steering effect. Normalization 
is made with respect to homogeneous flux. Solid circle: De- 
position flux. Open circle : Deposition flux after subtracting 
the enhancement due to the purely geometrical contribution 
(see the main text for details). 



III. RESULTS AND DISCUSSION 

As a preliminary investigation of the steering effect 
on thin film growth, the deposition flux distributions or 
deposition probabilities are examined for various depo- 
sition angles by MD. Deposition angle is measured from 
the normal to the (0,0,1) terrace to [-110] direction (y- 
axis). The positive deposition angle is for the deposition 
direction from the upper terrace to the lower one along 
y-axis as shown in Fig. 2(b), and the negative angle is 
for the opposite direction as shown in Fig. 2(a). The 
trajectories of the deposit atoms in Figs. 2(a) and (b) 
show the steering effect, where bending of trajectories of 
the incident atoms, most notably near steps, occurs due 
to the interaction between the deposit atom and sub- 
strate atoms. Figs. 2(c) to (g) show the deposition flux 
distributions normalized to homogeneous flux for various 
deposition angles. Depositon flux, shown with solid cir- 
cles in Figs. 2(c) to (g), increases near step , while that 
on terrace decreases compared with homogeneous flux?i^ 
It is important, however, to note that this enhanced flux 
near step is not soley due to the steering effect. The de- 
position flux distribution in Fig. 2(h) is for deposition 
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FIG. 3: Snapshot images of the 5ML Cu grown on Cu(l,l,17) 
at 240K. (a) Deposition with no steering effect, (b) Steered 
deposition at 70°. The size of the figures, (a) and (b), is 
800x150 ag. (c) Evolution of a step edge with increasing 
coverage. Successive curves show the development of step 
edges at Cu coverages below 5 ML with the increment of 1/3 
ML. The size of the figure (c) is 800x54 Oq. ao is the surface 
lattice constant of Cu(OOl), 2.55 A. 



with no steering effect considered, and still shows rela- 
tively high flux near steps. This is because there are only 
two adsorption sites available in 2.5 ag distance from each 
step edge along y-axis, while one adsorption site is avail- 
able in each 1.0 oq distance on terrace. The deposition 
flux after subtracting this purely geometrical contribu- 
tion is shown with open circles in Figs. 2(c) to (g) and 
shows the enhanced deposition flux near steps purely due 
to the steering effect. 

For deposition angles closer to the grazing angle (that 
is, anlges of larger magnitude), the deposition flux be- 
comes more inhomogeneous or more enhanced near steps, 
as can be seen by comparing Figs. 2(c) and (d) with Figs. 
2(e) and (f), respectively. As deposition angle becomes 
larger, so does the flight time of depositing atoms, dur- 
ing which their trajectories and in turn, the deposition 
fluxes are apt to be more disturbed by the inhomoge- 
neous substrate potential. It is also interesting to note 
the difference between the flux profiles at positive depo- 
sition angles (Figs. 2(d) and (f)) and those at negative 
deposition angles (Figs. 2(c) and (e)). In the negative 
angles, the deposition flux at the ascending step edge is 
larger than that at the descending step edge, and vice 
versa. This may be explained from the fact that at posi- 
tive deposition angles, the shadowing effect^" diminishes 
the probability for deposit atoms to sit on the sites next 
to the ascending step edge, while no such shadowing is 
expected for negative deposition angles. 

The effect of the steering-induced inhomogeneous de- 
position flux on thin film growth on a vicinal surface is 
studied by KMC utilizing MD for each deposition event. 
During the growth, the substrate temperature is set to 
240 K. Snapshots of a simulated system are shown in 
Figs. 3(a) and (b). Fig. 3(c) shows the evolution of a step 
as the coverage increases. We observe that the average 
position of step edge proceeds 8.5 oq for each monolayer 



(ML) deposition, indicating step fiow growth. However, 
the lateral roughness increases, and the coherence be- 
tween adjacent step edges develops to form 'finger'-like 
structures (Fig. 3(c)) as the coverage increases. Each 
finger shows ledge envelope along [100] and [0,-1,0] direc- 
tions, as observed for both experimental studies^ and the 
simulation results by Koponen et al..^ 

For a quantitative understanding of the growth insta- 
bility on a vicinal surface, the lateral roughness and the 
finger width taken as a measure of lateral coarseness are 
calculated. We define the lateral height, h{x), as the dis- 
tance from a position x at a pristine step edge to the 
growth front in the direction normal to the step edge 
(that is, in y-direction), and the lateral roughness as 

w{x) = \J < h{x)^— < h{x) >2 >. The lateral coarseness 
is calculated from the average separation between fingers 
within heights havg i 5ao, where havg is the average lat- 
eral height of each step. As a measure for the growth 
instability, we take the aspect ratio, lateral roughness to 
finger width (lateral coarseness). For an ideal step-flow 
growth or a stable growth, the aspect ratio should be 
very small. 

In Fig. 4(a), the lateral roughness increases monotoni- 
cally as coverage increases. At the maximum coverage of 
the present simualtion, 5 ML, the roughness is about 7ao 
indicating a very rough step edge. The roughness shows 
distinct dependence on the deposition angle. In the inset 
of Fig. 4(a), shown is the lateral roughness as a function 
of deposition angle after depositing 5 ML. The roughness 
is minimum at deposition angles at 0°. As the deposition 
angle becomes larger, so does the roughness. In addition 
to the deposition angle, the roughness depends also on 
the direction of deposition. When deposition is made fac- 
ing ascending step edge or at negative angle, the rough- 
ness of the fllm is small compared with that grown at the 
same magnitude of depostion angle, but in the opposite 
direction facing descending step edge. 

The development of lateral coarseness with increas- 
ing coverage was estimated by that of flnger width. In 
Fig. 4(b) and its inset, the finger width monotonically 
dereases as coverage increases, and also shows a definite 
dependence on both deposition angle and direction. The 
finger width shows maximum at —35° and decreases to 
minimum at +70°. The most notable thing is that the 
lateral roughness and coarseness have close correlation in 
their dependence on the deposition angle ; deposition at 
angles between —35° and 0° shows most stable step-flow 
growth with the minimum roughness and maximum fln- 
ger width or the minimum aspect ratio, while deposition 
at -1-70° shows the opposite behavior, the most unstable 
growth with the maximum roughness and the minimum 
flnger width, or the maximum aspect ratio. 

The aforementioned angular dependence of the growth 
instability should have originated from the dynamic ef- 
fect of deposition process, the steering effect, since all 
the kinetic variables are identical for each deposition at 
various angles. A direct result of steering effect is the 



4 



(0 
(A 
0) 

c 

JC 
D) 
3 
O 

75 

0) 

(3 



. — 0° (a) 

• +35° -ip^ 

-35° 

• — +70° 

70° 


NS ^s>^2 
,jr Ik 
y i 7.0 

/ ° 6.8 
/ to 

" / "* 
• / 5 6-6 


• 

• 

• . • 






-70 -35 §5 70 
Angle ( .) 



12 3 4 
Coverage (ML) 



■o 
c 




12 3 4 
Coverage (ML) 



FIG. 4: (a) Lateral roughness and (b) finger width (lateral 
coarseness) as function of coverage in the growth of Cu on 
Cu(l,l,17) at 240K. Refer main text for the definitions of lat- 
eral roughness and finger width. Inset: (a) lateral roughness 
and (b) finger width as a function of deposition angle after 
depositing 5 ML. The dotted lines in the figures (NS) and in- 
sets are the results of growth without considering the steering 
effect. 



inhomogeneous deposition flux. Hence, we investigate 
the correlation between depsoition flux distribution and 
growth stability: The atoms deposited near ascending 
step edge is expected to reproduce the step edge by di- 
rectly adhering to the sites near step edge, and should not 
be the main source of steering-induced growth instabihty. 
However, the atoms near the descending step edge would 
diffuse across terrace before reaching ascending step edge 
due to step Erlich Schwoebel barrier. During such terrace 
diffusion, the atoms redistribute themselves to feed and 
newly form laterally inhomogeneous structures, being a 
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FIG. 5: The aspect ratio of lateral roughness to finger width 
(solid circle) and the normalized deposition fiux averaged over 
three adsorption sites next to the descending step edge (open 
circle) are plotted as a function of deposition angle. Normal- 
ization is made with respect to homogeneous flux. 



source of meandering instability.* Indeed, we flnd the pre- 
dicted correlation between the growth instability and the 
enhanced deposition flux near descending step edge; In 
Fig. 5, the deposition flux averaged over the three sites ad- 
jacent to the descenting step edge is well matched with 
the aspect ratio for varying deposition angles. As the 
avearge deposition flux near the descending step edge is 
more enhanced, the mean travel length of deposit atoms 
to ascending step edge should become longer, and the 
growth becomes more unstable giving the larger aspect 
ratio. 

For possible origin of growth instability on a vicinal 
surface, two pictures have been proposed based on kinet- 
ics of adatoms; one attributes the instability to the step 
Erlich Schwoebel barrier effect (SESE)^:-S- and the other 
to the kink Erlich Schwoebel barrier effect (KESE).'^s^ 
SESE affects the motion and redistribution of deposited 
atoms on terrace, which should be directly dependent on 
the deposition flux destribution. KESE, however, gov- 
erns the motion of atoms along step edges, and is not 
directly affected by the initial deposition flux. The inti- 
mate correlation between deposition flux near descending 
step edge and growth instability shown in Fig. 5, indicates 
that the steering-induced deposition flux enhancement 
near descending step edge strengthens the role of SESE 
on growth instability. 

In the Figs. 4(a) and (b), the steered growth always 
show larger roughness and smaller coarseness regardless 
of the deposition angle than the growth neglecting the 
steering effect (dotted curves). That is, the steering effect 
always increases the growth instability. Such behavior is 
expected from the relatively small flux enhancement near 
descending step edge for steering- free depsoition as shown 
in Fig. 2, consistent with the aforementioned explanation. 
Although the steering effect is inevitable for vapor depo- 
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sition for thin film growth, the existence of a deposition 
angle producing the minimum growth instability (Fig. 5) 
suggests that the optimizaton of deposition angle should 
be a prerequisite for the most stable growth of thin films 
on a vicinal surface. 

IV. SUMMARY AND CONCLUSION 

KMC simulation in conjunction with MD simulation is 
performed to study the steering effect, in which the tra- 
jectory of each deposit atom is affected by interactions 



with substrate, on the growth of Cu on Cu(l,l,17). It is 
found that the steered deposition flux becomes inhomoge- 
neous and the inhomogeneity depends on the deposition 
angle and direction. The deposition flux enhancement 
near descending step edges is found to be the most criti- 
cal factor for the increase of growth instability due to the 
steering eflfect. The mechanism of such steering-induced 
increase of growth instability is discussed in details. In 
the present simulation, we also find a deposition angle 
producing minimum growth instability and show that the 
optimization of deposition angle should be a desirable for 
the most stable thin film growth on a vicinal surface. 
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